Exploiting geometric degrees of freedom in topological quantum computing 
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In a topological quantum computer, braids of non-Abelian anyons in a (2-|-l)-dimensional space- 
time form quantum gates, whose fault tolerance relies on the topological, rather than geometric, 
properties of the braids. Here we propose to create and exploit redundant geometric degrees of 
freedom to improve the theoretical accuracy of topological single- and two-qubit quantum gates. 
We demonstrate the power of the idea using explicit constructions in the Fibonacci model. We 
compare its efflciency with that of the Solovay-Kitaev algorithm and explain its connection to the 
leakage errors reduction in an earlier construction [Phys. Rev. A 78, 042325 (2008)]. 



I. INTRODUCTION 

Topological quantum computation is a rapidly devel- 
oping subject in recent years [l|, 0, H, [3|. In this novel 
scheme of quantum computation, information is stored 
in topological quantum states and intrinsically protected 
from local noises, and manipulation of quantum informa- 
tion is achieved by topological operations. A prototypical 
topological quantum computer is envisaged to be a sys- 
tem of exotic quasiparticles called non-Abelian anyons, 
which are believed to exist in various two-dimensional 
quantum systems H S S S i, 0, [H, [H [S Q M, 
[H, |T3]. A multiple of these anyons with fixed coordi- 
nates span a multi-dimensional Hilbert space, which can 
be used to construct qubits or encode quantum informa- 
tion [H, [l^ [2^ . Schemes have been propos ed to control 
and move anyons microscopically [2ll. l22l l23j . The world- 
lines of these anyons intertwine in (2-|-l)-dimcnsional 
space-time forming braids, which are quantum gates for 
topological quantum c omp utation. 

In earlier studies [l^, Ha, [2^ [l^l , researchers developed 
the method of brute-force search (and its variant) among 
braids within given braid length (measured by the num- 
ber of exchanges) to achieve a generic single-qubit quan- 
tum gate in the Fibonacci anyon model (which may be re- 
alizable in fractional quantum Hall systems [12, 2^]), and 
then constructed controUed-rotation gates from single- 
qubit gates. These works explicitly demonstrated the 
equivalence between a specific theoretical realization of 
topological quantum computer and a universal quantum 
computer model 0]. In general, a single-qubit quantum 
gate can be represented by a 2 x 2 unitary matrix 



G 



(1) 



where 6, a, (3 and 7 are real parameters. Apart from 
the overall phase factor e'", one needs three parameters 
6, (3 and 7 to specify the matrix. Within a given braid 
length, there are only a finite number of topological quan- 
tum gates, which form a discrete set in the U{2) space, 
thus generic gates can only be realized with a distribu- 
tion (wide on logarithmic scale) of error even in the ideal 
scenario (without technical or practical hindrance), due 



to the discrete nature of braid topology. This contrasts 
to many proposals of conventional quantum computa- 
tion, where quantum gates can be realized by contin- 
uously tuning physical parameters so generic quantum 
gates are expected to be realized with only a narrow dis- 
tribution of error (due to technical imperfections). On 
the other hand, the discreteness (thus error) in the re- 
alization of quantum gates with braids of finite length 
shares the same origin as the fault tolerance of topological 
quantum computation, as quantum states and quantum 
gates (braids) are topological and robust against local 
perturbations. This therefore poses an interesting ques- 
tion: How can we efficiently find the braid with finite 
length that approximates a desired quantum gate with 
error as small as possible? 

In a recent work, the authors proposed a novel con- 
struction of low-leakage topological quantum computa- 
tion based on the principle of error reduction by error 
introduction [26j . In topological quantum computation, 
the leakage errors in two-qubit gates one wants to mini- 
mize is often of topological origin, e.g., arising from the 
existence of noncomputational states. Nevertheless, one 
may find a class of equivalent braid constructions char- 
acterized by an additional geometric degree of freedom 
in braid segments, which can be said to correspond to 
a U{1) symmetry of the construction. In practice, how- 
ever, due to the discreteness of braids in the target space, 
such a symmetry is merely a pseudo symmetry. One finds 
that some of the constructions can have exponentially 
smaller errors than others ~ they are exactly what we 
want to find. The successful application of the princi- 
ple in the Fibonacci anyon model led to the discovery of 
an exchange braid (with a length of 99) that exchanges 
anyons between two different qubits, which can be used 
to construct generic controUed-rotation gates with leak- 
age error as small as 10~^. However, the idea can not be 
directly applied to construct single-qubit gates because 
apparently there is no such geometric freedom. 

In this paper, we generalize the idea that errors in 
topological gate construction can be reduced by the intro- 
duction of an additional geometric redundancy (or sym- 
metry) ^ and show that it can also be applied to the 
construction of generic single-qubit gates with unprece- 
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dented efficiency and accuracy in theory. We demon- 
strate this idea exphcitly in the Fibonacci model, though 
it is apphcable in generic models that support universal 
topological quantum computation. By introducing new 
degrees of freedom with unitary similarity transforma- 
tion, we show a generic single-qubit gate can be approx- 
imated to a distance [29^ of the order 10~^° by a braid 
of length ~300, which is more efficient than a direct ap- 
plication of the Solovay-Kitaev algorithm [2^. We also 
discuss the significant reduction of leakage error by er- 
ror introduction in a parallel construction of two-qubit 
controUed-gates (see also Ref. [13, H^) to demonstrate 
the power and generality of the idea of topological error 
reduction by exploiting redundant geometric degrees of 
freedom. 



II. SINGLE-QUBIT GATE CONSTRUCTION 

Let us first discuss the high-accuracy construction of 
a generic single-qubit gate G represented by Eq. To 
create additional degrees of freedom needed for error re- 
duction, one is tempted to decompose the target gate as 
G = G1G2, where 



Gi,2 = e'"i'^ 



1 - 62_2 e-"^''^ 
-6i,2e-^^i'^ 



1 



/,2 pi/3i,2 
"1.2 



(2) 



are unitary matrices, or quantum gates. Unfortunately, 
the degrees of freedom of the two gates are dependent, 
as can be immediately seen from G\ — GG2 , i.e., one 
of the gate (e.g., Gi) is completely fixed by the other 
gate (e.g., G2) up to an unimportant global phase factor. 
Therefore, one cannot use the new degrees of freedom for 
gate optimization. In fact, this decomposition of G is the 
mathematical structure of the bidirectional search [2^ . 
which facilitate the search for twice longer braids without 
a significant increase in CPU power. 

In fact, to create separable degrees of freedom with two 
gates, we can write G = GiG2G^, a unitary similarity 
transformation, which creates geometric redundancy. It 
is easier to visualize the transformation in terms of rota- 
tion in three dimensions, thanks to the homomorphism 
between the groups 5*0(3) and SU{2). This means that 
a rotation around an arbitrary axis / by an angle on 
a Bloch sphere can be carried out by first rotating / to 
another direction I', then rotating around I' by an an- 
gle 9, and finally rotating /' back to /. The geometric 
interpretation clearly indicates that the freedom in the 
choice of can be exploited to optimize the single-qubit 
gate G, because these apparently equivalent realizations 
use different braids to approximate the target gate with 
entirely different accuracy. 

We can use a phase gate P (a diagonal matrix) to 
illustrate the determination of Gi and G2 without loss 
of generality. This is because, according to the spectral 
theorem for normal matrices, any unitary matrix G can 
be unitarily diagonalized as G = S^PS, where 5 is a 



unitary matrix. We can then contract S with Gi, so 
that P = G1G2G+, where Gi = ^Gi. 
For concreteness, let us assume 



P = 



■-'^ 




(3) 



which is a rotation around the z axis by an angle (3. 
The parameters 61^2, 01.2 and 71^2 of Gi and G2 that 
decompose P must, therefore, satisfy 



bi 



(1 - &^)^/2cos/32 cos/3, 

&2 



2 sin^ /? + 2(1 - 6^)1/2 sin /32 sin /3 
/3i + 71 = 72 + (fc + l/2)7r. 



(4) 
(5) 
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where the integer k is even for positive sin/3 or odd for 
negative sin/3 (we exclude the case sin/3 = when the 
corresponding gate is proportional to the identity) . From 
Eq. (HI we can see that G2 is a rotation of the same angle 
as that of P but around a new axis, which is related to 
the original rotation axis of P by Eqs. ([5]) and Never- 
theless, when G2 is fixed, Gi is only partially determined 
by G2 and P and still has a degree of freedom (between 
f3i and 71). In other words, the similarity transformation 
G1G2G1 has an SU{2) symmetry, from which we have 
three free parameters to choose. Gi, the rotation of z- 
axis to a new axis, has a U{1) symmetry (i.e., one free 
parameter), while G2, fixing the direction of the new axis, 
has a symmetry of SU{2)/U{1) ^ S"^ (i.e., two free pa- 
rameters). Hence, we can successfully separate the three 
degrees of freedom into two parts in Gi and G2, which al- 
low us to efficiently search for high-accuracy single-qubit 
gates. 

As an explicit demonstration of the algorithm, we con- 
struct a phase gate 
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in the Fibonacci anyon model (please refer to Refs. [26|, 
for details of this model), where there are two types 
of anyons with topological charges (vacuum) and 1 (Fi- 
bonacci anyon) satisfying a nontrivial fusion rule 1x1 = 
0-1-1. We use two pairs of Fibonacci anyons with to- 
tal charge to encode one bit of quantum information. 
The basis states are chosen as |0) — |((ll)o(ll)o)o) and 
|1) = |((ll)i(ll)i)o), where the subscripts specify the fu- 
sion results (or total topological charges) of the anyons 
in the preceding brackets. Therefore, four-strand braids 
can be generated by the elementary braids with repre- 
sentation 



0-1 = 0-3 = 
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FIG. 1: (a) Mapping two qubits (ai-a^ and as-ag) into one 
qubit consisting of four composite anyons A1-A4. (b) Braid 
P2 of the composite anyons Ai to A4. As explained in the 
text, we only move the composite anyon A2 to braid with 
the composite anyons A3 and A4 and return A2 back to the 
original position at the end of the braid. 



and their inverses, where r — (\/5 — We find a 
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verify G1G2G1 , with 280 interchanges [3^J, approxr 
mates Pi with a distance '--^ 4 x 10^^". In general, such a 
precision can be achieved by a braid of length ~300 for 
a generic single-qubit gate with the algorithm specified 
above. A braid of similar accuracy is expected to exist 
at a length of as short as 150 [2^, but the search for it 
is exponentially harder. 



III. TWO-QUBIT GATE CONSTRUCTION 

The single-qubit construction scheme echos the earlier 
low-leakage two-qubit construction scheme [2^ , in which 
one exchanges a two-anyon composite in the control qubit 
with the neighboring anyon in the target qubit and per- 
forms single-qubit operations on the new target qubit, 
which translate into controlled rotations in the original 
two-qubit system, before one exchanges the anyons back 
to their original locations. There is, however, a notable 
difference: the SU (2) symmetry in the single-qubit con- 
struction is broken in the two-qubit construction due to 
leakage error. As discussed in the introduction, only a 
U{1) symmetry exists in the two-qubit construction when 
we enforce leakage errors to be negligible. This suggests 
that the requirement of zero (or ultralow) leakage errors 
in the two-qubit construction eats two degrees of free- 
dom. It becomes clear in an alternative high-accuracy 
(^10^^*^) implementation presented in the following, in 
which the construction of the two-qubit gates are based 
on a mapping from two qubits to one qubit in the four- 
anyon encoding scheme (see also [13, &^), which also 
demonstrates the generality of the idea of topological er- 
ror reduction by redundant geometric degree of freedom. 

Explicitly, the mapping scheme is the following. For 
clarity, we label the anyons in the target qubit 01-04 and 
those in the control qubit 05-03 as in Fig. [Ija). We can 
treat the two pairs of anyons in each qubit as two com- 



posite anyons, which have a total topological charge 0. 
Then we have a mapping from two qubits of Fibonacci 
anyons to one qubit of composite anyons, which we label 
A1-A4 as in Fig. [T](a). The computational basis states 
are chosen as 

|00) = |((ll)o(ll)o)o((ll)o(ll)o)o) = |((00)o(66)o)o), 
101) = |((ll)o(ll)o)o((ll)i(ll)i)o) = |((00)o(ii)o)o), 

|io) ^ |((ii)i(ii)i)o((ii)o(ii)o)o> = |((ii)„(66)o)o), 
111) = |((ii)i(ii)i)o((ii)i(ii)i)o) = |((ii)o(ii)o)o), 

(10) 

where and 1 denote the topological charge of compos- 
ite anyons. In fact, the composite qubit is not a qubit 
in the normal encoding scheme, because the composite 
anyons Ai can have charge 0. Each pair of composite 
anyons (e.g., Ai and A2) always have total charge 0, un- 
less leakage error occurs so each of the original qubits 
has total charge 1. Therefore, the class of braids that 
we look for to manipulate the the composite qubit as in 
Fig. [T](b) without introducing leakage errors are the ones 
that realize phase gates, for example. 
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(11) 



We see immediately that two degrees of freedom in SU (2) 
disappear due to the requirement of zero leakage errors 
and we are left with a U{1) symmetry only. In prac- 
tice, we restrict ourselves to move the composite anyon 
A2 to braid with the composite anyons A3 and A4 and 
return A2 back to the original position at the end of the 
braid; this is known as a weave [l^]. When the compos- 
ite qubit is in the state |((ll)o(ll)o)o), this braid will 
introduce a phase factor e^^°'^~^^^ to the two-qubit sys- 
tem [e.g., a phase factor -1 for the braid approximating 
Pi in Eq. ([7])]. While if the composite qubit is originally 
in the other computational states in Eq. PU)) . either the 
topological charge of the composite anyon A2 or the topo- 
logical charges of the composite anyons ^3 and A4 are 
0. Since the braid between an anyon with topological 
charge and another anyon with topological charge ei- 
ther or 1 does not change the state of the system, the 
braid will bring only a trivial phase factor 1 to the system. 
Thus a braid approximating the single-qubit phase gate 
P2 in Eq. pT|) corresponds to a controUed-phase gate, 
e.g., a controUed-Z gate for the braid approximating Pi 
in Eq. dZD. 

A scheme to construct an arbitrary controUed-rotation 
gate, parallel to the single-qubit gate construction, is il- 
lustrated in Fig. [2] (see also Fig. 3 in Ref. [s^ for an 
SU{2)5 construction). In particular, we need to ap- 
ply a single-qubit gate G3 on the target qubit after the 
controUed-phase gate and its inverse G3"^ = G+ before 
the controUed-phase gate. This is in the same spirit as 
the similarity transformation in the single-qubit case, ex- 
cept that the controUed-phase gate is defined on the com- 
posite qubit, not the target qubit. For completeness, we 
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FIG. 2: A scheme to realize a generic controlled-rotation gate 
in the Fibonacci anyon model. We consider two qubits com- 
posed of ai-as- P2, which realizes a two-qubit controUed- 
phase gate, is a braid acting on the effective single qubit 
formed by composite anyons A1-A4. G3 and P3 are braids 
of single-qubit gates that modify the controUed-phase gate to 
a generic controlled-rotation gate. Note that we choose the 
convention of the time direction from left to right. 



need to introduce another single-qubit phase gate 
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to adjust the phase of the control qubit such that the 
resulting gate is 
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where R is related to G3 by 



^+2&)G3 



i(a2-/32)/2 Q 

Q gi(a2-/52)/2 



G3+. (14) 



Given a target i?, we should search for a braid realizing 
the corresponding G3. As R can be diagonalized by a 
similarity transformation SRS~^, the constraint on G3 
is that SG3 should a single-qubit phase gate with an 
arbitrary phase, which again allows a redundant U{1) 
degree of freedom. 

In fact, we intentionally designed our presentation 
so that the braid realization of Pi in Eq. ([7]) is the 
same for P2 in order to construct a controlled-NOT 
(CNOT) gate. Correspondingly, we find the sequence 
for one instance of G3 as a 



—4 —4 2 2 — 



-4.^2 
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a^O^ CT2Cr3 U2"3 



o^'^a\a^^a\a\a2'^a\a\a\a\. The total braid for the 
CNOT gate, with an error of 5 x 10~^", contains 280 
exchanges of double braids and 208 exchanges of single 
braids. Note that P3 is trivial for the CNOT gate. 
We would like to point out that this construction is 
conceptually interesting but technically less efficient, 
because an exchange of two double braids is, in fact, 
four exchanges of single braids. The new single-qubit 
construction combined with the two-qubit construction 
in the earlier proposal (26j , which also involves exchanges 
of a single braid and a double braid, can achieve a similar 
error within 1000 exchanges of single braids. 
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IV. CONCLUSION AND DISCUSSION 

In conclusion, we proposed the idea of exploiting re- 
dundant geometric degrees of freedom in topological 
quantum computation to reduce the topological errors 
due to discreteness of gates realized by finite-length 
braids. This is possible because we can separate the re- 
dundant degrees of freedom into (partially) independent 
parts, which allows topological quantum gate construc- 
tion to be more efficient. We also established the intrigu- 
ing connection between the sacrifice of two such degrees 
of freedom and the minimization of two-qubit leakage er- 
rors. 

We can understand the error reduction from a different 
angle. In the three-dimensional space of unitary matri- 
ces, a target gate is just a zero-dimensional point. The 
introduction of geometric redundancies transforms the 
target into a one- or higher-dimensional object, thereby 
allows an efficient deeper search. The algorithm is prac- 
tically useful as computational errors can be reduced ex- 
ponentially at all length scales by the introduction of 
redundant degrees of freedom as shown, e.g., in Fig. 4 
of Ref. [2^ . The increase in braid length by a factor of 
roughly three is thus well conpcnsated by the exponential 
suppression in error. 

Finding optimal braids belongs to the generic ques- 
tion of approximating an arbitrary unitary operation by 
a set of discrete gates (or matrices) relevant to, e.g., con- 
structing quantum circuits in generic quantum compu- 
tation, for which a remarkable rate of convergence can 
be achieved by the Solovay-Kitaev algorithm [3^. The 
Solovay-Kitaev algorithm is based on the principle of 
error cancellation by a group commutator structure of 
ABA^^B^^ factor (given an initial e-net). It has been 
implemented in the context of topological quantum com- 
putation by Hormozi et al, who achieved a gate with a 
distance ~ 4.2 x 10~^ to iX with a braid of length 220 
in one iteration [25j . Therefore, the algorithm presented 
here can achieve comparable accuracy to that from apply- 
ing one iteration of the Solovay-Kitaev algorithm, albeit 
with braids that are about 40% shorter than those ob- 
tained from the Solovay-Kitaev algorithm. An iterable 
modification of the algorithm, as well as its performance 
comparison with the Solovay-Kitaev algorithm, is pre- 
sented elsewhere [33. 
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